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\Z!^] Abstract. The conjectures of Deligne, Beihnson, and Bloch-Kato as- 

'^1^ . serf that there should be relations between the arithmetic of algebro- 

I — ^ ' geometric objects and the special values of their L-functions. We make a 

numerical study for symmetric power L-functions of elliptic curves, ob- 
taining data about the validity of their functional equations, frequency 
of vanishing of central values, and divisibility of Bloch-Kato quotients. 

H 

• ■ 1 Introduction and Motivation 

C^ ' There are many conjectures that relate special values of L-functions to the 

arithmetic of algebro-geometric objects. The celebrated result C(2) — 7r^/6 of 
Euler [2(J[ §XV] can be reinterpreted as such, but Dirichlet's class number for- 
mula |15[ §5] is better seen to be the primordial example. Modern examples run 
^ ■ the gamut, from conjectures of Stark [lU] on Artin L-functions and class field 

l/^ I theory, to that of Birch and Swinnerton-Dyer |2] for elliptic curves, to those of 

On ■ Beilinson |1I32| related to i^-theory, with a passel of others we do not mention. 

^P I For maximal generality the language of motives is usually used (see |21 §1-4]). 

One key consideration is where the special value is taken. The L-function 
vQ I can only vanish inside the critical strip or at trivial zeros; indeed, central values 

f^ ■ (at the center of symmetry of the functional equation) are the most interesting 

ones that can vanish, and the order of vanishing is likely related to the rank of a 
geometric object (note that orders of trivial zeros can be similarly interpreted) . 
We have chosen to explore a specific family of examples, namely symmetric 
power L-functions for rational elliptic curves. The impetus for this work was 
^ I largely a theoretical result ^^1 of the first author, whose computation of Euler 

factors in the difficult case of additive primes greatly reduced the amount of 
. ^ hassle needed to do large-scale computations. Previous theoretical work includes 

j^ ■ that of Coates and Schmidt [7] on the symmetric square and Buhler, Schoen, and 

Top on the symmetric cube W ; this second paper also contains a lot of compu- 
tational evidence, while Watkins has provided much data 44 in the symmetric 
square case. In some cases the L-functions we use are not known to possess the 
properties that would be required to justify that our computations produce num- 
bers of any validity whatsoever — in these cases, the "numerical coincidence" 
in our computations can be seen as evidence for the relevant conjectures. 

* Supported during parts of this research by EPSRC grants GR/N09176/01 and 
GR/T00658/01, the Isaac Newton Institute, the CNRS and the Institut Henri 
Poincare, and the MAGMA Computer Algebra Group at the University of Sydney. 
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2 ly- functions 

We define the symmetric power L-functions of an elliptic curve E/Cl via comput- 
ing an Euler factor at every prime p. This Euler factor is computed by a process 
that essentially just takes the symmetric power representation of the standard 
2-dimensional Galois representation associated to E, and thus our method is a 
generalisation of that used by Coates and Schmidt [7' for the symmetric square/ 
following the original description of Serre [36| . We briefly review the theoretical 
framework, and give explicit formulae for the Euler factors in a later section. 

For every prime p choose an auxiliary prime I ^ 2,p and fix an embedding 
of Qi into C. Let Et denote the t-torsion of E, and Ti{E) — liraEin be the 

/-adic Tate module of E (we fix a basis). The module Vi{E) — Ti{E) (g)z, Q; 
has dimension 2 over Q; and has a natural action of Gal(Qp/Qp) [indeed one of 
of Gal(Q/Q)], and from this we get a representation /?/ : Gal(Q /Qp) — > Aut(V/). 
We write Hl[E) = Homq, (V;(i?), Q;), and take the mth symmetric power of 
the contragredient of pi , getting 

pr : Gal(Qp/Qp) -> Aut(Sym" (///(£;))) C GL™+i(C). 

We write Dp = Gal(Qp/Qp), let Ip be the inertia group of this extension, and 
let Frobp be the element of Dp/Ip = Gal(Fp/Fp) given by x ^ x^ . With all of 
this, we have 

L(Sym"i?,s) = []det [ld„,+i - pT{Yvoh-')p-^ | (Sym™(i7i (£;)))'' 
p 

For brevity, we write Lm{E,s) = i(Sym™£', s), and denote the factors on the 
right side by Um{p;s). As mentioned by Coates and Schmidt [3 p. 106], it 
can be shown that Um{p',s) is independent of our choices. The analytic theory 
and conjectures concerning these symmetric power L-functions are described 
in inni- In particular, the above Euler product converges in a half-plane, and is 
conjectured to have a meromorphic continuation to the whole complex plane. 

We also need the conductor A,„ of this symmetric power representation. 
We have A,„ = JlpP^"^^^ where /„(p) = e,„(/p) -I- S^ip). Here em(/p) is the 

codimension of (Sym"(i//(£;)))'^'' in Sym'"(iJ/(£;)); we shall see that it can 
be computed via a character-theoretic argument. The wild conductor Smip) is 
unless p = 2, 3, when it can be computed as in ^^ §2.1] or the appendix of [7]- 

^ Note that Buhler, Schoen, and Top |^ phrase their definition of Euler factors differ- 
ently, as they emphasise that conjecturally the L-function is related to a motive or 
higher-dimensional variety; however, their definition is really the same as ours. 



2.1 Critical values 

The work of Deligne ^J Prop. 7.7S\ tells us when and where to expect critical 
values; these are a subset of the more-general special values, and are the easiest 
to consider.^ When m = 2v with v odd there is a critical value LmiE,v + 1) 
at the edge of the critical strip, and when m = 2u — 1 is odd there is a critical 
central value Lra{E,u). We let i?+, i7_ be the real/imaginary periods of E for 
m = 1, 2 (mod 4), and vice- versa for to = 3 (mod 4). In the respective cases of 
TO even/odd we expect rationality (likely with small denominator) of either 

L.^{E,v + l)f 27rN y^^+''^/^ _ L,„(ii;,«)(2^iV)«("-i)/2 



(27r)"+l \f2+[2^J ^u(u+l)/2^u(u-l)/2 

When m is odd, the order of L„i {E, s) at s — u should equal the rank of an associ- 
ated geometric object. The Bloch-Kato conjecture ^ relates the quotients in |^ 
to iJ^-groups, Tamagawa numbers, and generalised Shafarevich-Tate groups.'^ 

3 Computation of Euler factors and local conductors 

We first consider multiplicative and potentially multiplicative reduction for a 
given prime p; these cases can easily be detected since Vp{JE), the valuation of 
the j-invariant, is negative, with the reduction being potentially multiplicative 
when p\c4. When E has multiplicative reduction, the filtration of [HJ §8] implies 
the local tame conductor £„ is m and 5m(2) = S„i{3) = for all m. The Euler 
factor is U.m{p', s) = (1 — a^/p^)'-^, where Up = ±1 is the trace of Frobenius. In 
the case of potentially multiplicative reduction, for m odd we have e™ == to, + 1, 
and so Um{p] s) = 1, while with m even, we have that em = in and compute that 
Umip; s) = (1 - l/p')~^. The wild conductor at p = 2 is 5„(2) = ^Si{2) for 
odd m and is zero for even to,, while (5m (3) = for all to,. 

3.1 Good and additive reduction — tame conductors 

Let E have good or potentially good reduction at a prime p, and choose an 
auxiliary prime I ^ 2,p. The inertia group Ip acts on Vi{E) by a finite quotient 
in this case. Let Gp = Gal(Qp(i?/)/Qp) and ^p be the inertia group of this 
extension.'' The work of Serre 37 lists the possibilities for <Pp. It can be a cyclic 
group Cd with d = 1,2, 3, 4, 6; additionally, when p = 2 it can be Qg or SL2{F3), 
and when p = 3 it can be C3 x: C4. For each group there is a unique faithful 
2-dimensional representation !Z',f of determinant 1 over C, which determines p;. 

^ Critical values conjecturally only depend on periods (which are local objects), while 
the more-general special values can also depend on (global) regulators from K-theoiy. 

^ See |18l §7] for an explicit example; note his imaginary period is twice that of our 
normalisation (and the formula is out by a power-of-2 in any case), and the conductor 
enters the formula in a different place (this doesn't matter for semistable curves). 

* The group <Pp is independent of the choice of I (see [371 p. 312]), while only whether 
Gp is abelian matters, and this independence follows as in Lemmata 1.4 & 1.5]. 



Our result now only depends on <?; for a representation ^ we have the trace 
relation (which is related to Chebyshev polynomials of the second kind) 

tr(Sym'"if ) ^ ^ ( ™ jtr(<^')"-2fe(_detl^)^ (2) 

and from taking the inner product of tr(Sym'"tf'<f ) with the trivial character we 
find the dimension of the ^- fixed subspace of Sym"^(^Hl{E)) , which we denote 
by Pm (^) ■ Upon carrying out this calculation, we obtain Table Q] which lists 
values for /3m{'^), from which we get the tame conductor em(^) = m+1 — /3m(<?). 
The wild conductors 5m{p) are for p > 5, and for p = 2,3 are described below. 

Table 1. Values of f3mi^) for various inertia groups; here fh is m modulo 12. 



m 


C2 


Cs 


Ci 


Ce, 


Qs 


C-i » C4 


SL2(F3) 





m + 1 


(m + 3)/3 


(m + 2)/2 


(m + 3)/3 


(m + 4)/4 


(m + 6)/6 


(m + 12)/12 


1 





(m-l)/3 

















2 


m + 1 


(m + l)/3 


m/2 


(m+l)/3 


(m-2)/4 


(m-2)/6 


(m-2)/12 


3 





(m + 3)/3 

















4 


m + 1 


(m-l)/3 


(m + 2)/2 


(m-l)/3 


(m + 4)/4 


(m + 2)/6 


(m-4)/12 


5 





(m + l)/3 

















6 


m + 1 


(m + 3)/3 


m/2 


(m + 3)/3 


(m-2)/4 


m/6 


(m + 6)/12 


7 





(m-l)/3 

















8 


m + 1 


(m + l)/3 


(m + 2)/2 


(m + l)/3 


(m + 4)/4 


(m + 4)/6 


(m + 4)/12 


9 





(m + 3)/3 

















10 


m + 1 


(m-l)/3 


m/2 


(m-l)/3 


(m-2)/4 


(m-4)/6 


(m-10)/12 


11 





(m + l)/3 


















3.2 Good and additive reduction — Euler factors for p > 5 

When p > 5, a result of Serre |2Z1 tells us that the inertia group is 'P = Cd where 
d = 12/ gcd[l2,Vp{AE)) ■ Note that this gives d = 1 when p is a prime of good 
reduction, which we naturally include in the results of this part. We summarise 
the results of Martin's work ^^ concerning the Euler factors. Note that the 
result of D§,browski QJl Lemma 1.2.3] appears to be erroneous. 

There are two different cases for the behaviour of the Euler factor, depend- 
ing on whether the decomposition group Gp = Gal(Qj,(£'i)/Qp) is abelian. 
From 123 Prop. 2.2] or ^A, Th. 2.1], we get that this decomposition group is 
abelian precisely when p = 1 (mod d) . When Gp is nonabelian we have 



Um{p;s) - (1 - {~pr^'/pT''-a + {-pr'vf) 



m/2 , 



(3) 



vG„ 



where Am + i?m = Pm and A,„ is the dimension of (Sym™(iJ;^ (£'))) " . Using 
Gp/^p = C2 and det{\I^,p{x)) = —1 for x e Gp\<?p, more character calculations 
tell us this dimension is (/3m + l)/2 when /3m is odd and is /3m/2 when Pm is even. 
This also holds for the non-cyclic 'P when p = 2,3, for which Gp is automatically 
nonabelian. When (p — C3 and m is odd, we have Um{p\ s) = (1 +p'^ lp^'')~^"'/'^. 



When Gp is abelian, we need to compute a Frobenius eigenvalue Up (whose 
existence follows from 38, p. 499]). In the case of good reduction, this comes from 

counting points mod p on the elliptic curve; we have ctp = (ap/2) ± i\ p — <ii,/^ 

where p -\- 1 — ap is the number of (projective) points on E modulo p. And 
when <P = C2 we count points on the pth quadratic twist of E. In general, we 
need to re-scale the coefficients of our curve by some power of p that depends 
on the valuations Vp of the coefficients. Since p > 5, we can write our curve 
as y'^ = x^ + Ax + B, and then re-scale by a factor t = pmin(«p(^)/2,t;p(B)/3) ^.^ 
get a new curve _E* : y^ == a;"^ + Ax/t'^ + B/t^, possibly defined over some larger 
field. Because of our choice of t, at least one of A/t'^ and B/t^ will have Vp equal 
to 0. The reduction E* modulo some (fractional) power of p is then well-defined 
and non-singular, and we get ap from counting points on E*; it turns out that 
choices of roots of unity will not matter when we take various symmetric powers. 
Returning back to Um{p] s), we get that when Clp{Ei) / Qp is abelian this Euler 
factor is 

Um{p;s)= n (i-«r^«;//)-'. (4) 



0<i<71 

d\(2i-r. 



3.3 Considerations when p = 3 

Next we consider good and additive reduction for p = 3. We first determine the 
inertia group, using the 3-valuation of the conductor as our main guide. In the 
case that V3{N) = we have good reduction, while when V3{N) = 2 and V3{A) is 
even we have "P — C2- Since G3 is abelian here, the Euler factor is given by Q, 
while the wild conductor is and tame conductor is obtained from Table Q 
When V3{N) = 2 and V3{A) is odd we have that <P = C4 and G3 is nonabelian. 
The wild conductor 5m (3) is 0, and the Euler factor is given by Q. 

When V3{N) = 4 we get # = C3 or Ce, the former case when 4|-y3(Z\). For 
these inertia groups, the question of whether G3 is abelian can be resolved as 
follows (see ^^ Th. 2.4]). Let 64 and ce be the invariants of the minimal twist 
of E at 3. In the case that C4 = 9 (mod 27), we have that G3 is abelian when 
C6 = ±108 (mod 243) while if 3^]c4 then G3 is abelian when £4 = 27 (mod 81). 
In the abelian case we have as = Ci2v3 up to sixth roots, which is sufficient. 
The Euler factor is then given by either (jJJ or Q, the tame conductor can be 
obtained from Table Q] and the wild conductor (computed as in the appendix 
of 13) from Table El When -i;3(iV) = 3, 5 we have that <? = C3 xi C4. The Euler 
factor is given by lj3l and the wild conductor can be obtained from Table|21 with 
the first C3 XI C4 corresponding to V3{N) — 3, and the second to V3{N) = 5. 

3.4 Considerations when p — 2 

Finally we consider p == 2, where first we determine the inertia group. Let M 
be the conductor of the minimal twist F of £' at 2, recalling 44, § 2.1] that in 
general we need to check four curves to determine this twist. Table0]then gives 
the inertia group. The appendix of [J] omits a few of these cases; see 03]. When 



<? = Ci , 6*2 we can always determine ap via counting points modulo p on E or a 
quadratic twist, and G2 is always abelian. The Euler factor is then as in Q). For 
<l^ — CsjCq the group G2 is always nonabelian, and the Euler factor is as in 0. 
For the case of <? = C4 and p — 2, the question of whether G2 is abelian comes 
down |44[ Th. 2.3] to whether the C4 invariant of F is 32 or 96 modulo 128, it 
being abelian in the latter case, where we have a2 = CsV^ up to fourth roots. 
The Euler factors for this and the two cases of noncyclic ^ are obtained from ||3J) 
or Q, while the wild conductors (5m(2) are given in Table[21 with the appropriate 
line being determinable from the conductor of the first symmetric power. 



Table 2. Values for (5„(2). 


#2 


m = 1 


formula 


C2,C6 


2 


em(C2) 


C2,Ce 


4 


2e™(C2) 


Ca 


6 


26,„(C4)+e™(C2) 


Qs 


3 


em(Q8) + iem(C2) 


Qs 


4 


imiQs) +e,n{C2) 


Qs 


6 


E™(Qs) + f™(C4) + em(C2) 


SL2(F3) 


1 


|e™(08) + |em(C2) 


SL2(F3) 


2 


lemiQa) + lemiCi) 


SL2(F3) 


4 


lemiQs) + |em(C2) 


SL2(F3) 


5 


fem(Q8) + |em(C2) 



Table 3 


Values for (5^(3). 


<P3 


m = 1 


formula 


C-j,, Cfi 

C3 XC4 
C3 XC4 


2 

1 

3 


2^m{Cz) 

|e™(C73) 



Table 4. Values of '?2• 


«2(Af) 


$2 



2 

3,7 
5 

8 


Ci if t;2(Af) =0elseC2 
C3 if «2(iV) =2elseC6 

S'L2(F3) 

gs 

Qs if 2''|C6(F) elseC4 



3.5 The case of complex multiplication 

When E has complex multiplication by an order of some imaginary quadratic 
field K, the situation simplifies since we have L{E, s) = L{tp, s — 1/2) for some^ 
Hecke Grossencharacter ip. For the symmetric powers we have the factorisation 



L(Sym™£:, s) 



m/2 



m — 2i 



m/2), 



(5) 



i=0 



where V" is the ^-function when A\m, and when 2||m it is L{9k,s) for the 
quadratic character 9k of the field K. Note that the local conductors and Euler 
factors for each L^ij)^ , s) can be computed iteratively from jSJl since this informa- 
tion is known for the left side from the previous subsections. This factorisation 
reduces the computational complexity significantly, as the individual conductors 
will be smaller than their product; however, since there are more theoretical re- 
sults in this case, the data obtained will often lack novelty. The factorisation ((SJ 
also implies that L2u-i{E, s) should vanish to high degree at s = u, since each 
term has about a 50% chance of having odd functional equation. We found some 
examples where L{ijj^,s), L{ip^,s), or L{ili' ,s) has a double zero at the central 
point, but we know of no such triple zeros. 

^ This is defined on ideals coprime to the conductor by il^{z) — xi\z\)(z/\z\) where z is 
the primary generator of the ideal and x is generally a quadratic Dirichlet character, 
but possibly cubic or sextic ii K = Q(\/— 3), or quartic ii K = Q(V~l)- When 
taking powers, we take X"* to be the primitive Dirichlet character which induces X'' ■ 



4 Global considerations and computational techniques 

We now give our method for computing special values of the symmetric power 
L- functions defined above. To do this, we complete the L- function with a r'-factor 
corresponding to the prime at infinity, and then use the (conjectural) functional 
equation in conjunction with the method of Lavrik |29) to write the special 
value as a "rapidly-converging" series whose summands involve inverse Mellin 
transforms related to the .T- factor. First we digress on poles of our L-functions. 

4.1 Poles of i- functions 

It is conjectured that Lm{E, s) has an entire continuation, except when 4|?Ti and 
E has complex multiplication (CM) there is a pole at s = 1 + m/2, which is the 
edge of the critical strip. ^ We give an explanation of this expectation from the 
standpoint of analytic number theory; it is likely that a different argument could 
be given via representation theory. We write each Euler factor as [/,„(p;s) = 

(l - brnip)/?" H ) and as s ^ 1 + m/2 we have log L,„(s) -- J^p bm{p)/p^- 

We will now compute that the conjectural Sato- Tate distribution ^41^ implies that 
the average value of bm{p) is 0, while for CM curves the Hecke distribution [25 
will yield an average value for bm{p) of p"'/^ when 4|to. 

Similar to Q, for a good prime p we have b^ip) = Y^T=o ("TO°^™~^*(~P)*- 
The Sato- Tate and Hecke distributions imply that the average values of the fcth 
power of Op are given by 

. , C (cos 6')'=^ (sin 6*) 2^6* . .ClcosOfde 

We have (aj;) — for k odd; for even k the Wallis formula 021 implies 
(COS0) (sm^) de - ^ (fc + 2)!! " (fc + 2)!!' 

(fc+2)!! 



SO that (Cp) is (2^)^ (fc+2')!!' • ^"^ induction exercise shows that this implies 
{bm{p)) = when E does not have CM. We also have {al)cM = {'^^f'^T^ 
for even k, and again an inductive calculation shows that {bm{p))cM = p'"'^ 
when A\m and is zero otherwise. This behaviour immediately implies the afore- 
mentioned conjecture about the poles of Lm{E, s) at s = 1 + m/2. 

4.2 Global considerations 

Let A„XE,s) = C^7„j(s)i,„(S,s), where C^ = A^„j/(27r)"+i for m odd and is 
twice this for m even. For m odd we write m = 2u — 1 and for m even we write 
m = 2v; then from J14l §5.3] we have respectively either 

u~l v—1 

lm{s)=\{r{s-l) or -lra{s)=r{s/2~[v/2\)\{r{s~i). 
1=0 i=0 



The case of m = 4 follows as a corollary of work of Kim |2(il Corollary 7.3.4]. 



When 4|m and E has CM, we multiply 7m(s) by (s — v){s — v — 1). We ex- 
pect Am{E,s) to have an entire continuation and satisfy a functional equation 
AmiE,s) = WmAm{E,m + 1 — s) for some Wm = ±1. The works of Kim and 
Shahidi jJT] establish parts of this conjecture J We can find Wm via experiment 
as described in Section 14.41 but we can also try to determine w„i theoretically. 

4.3 Digression on local root numbers 

The sign Wm can theoretically be determined via local computations as in |12| . 
but this is non-trivial to implement algorithmically, especially whenp = 2, 3. We 



expect to have a factorisation Wr, 



n. 



w. 



i{p) where the product is over bad 



primes p including infinity. For m even, the very general work of Saito |35| can 
then be used to show^ that Wm = -Hi, so we assume that m is odd. From |14[ §5.3] 
we have u'm(oo) = —{-j^)', combined with the relation Wmip) = wi{p)"^ for 
primes p of multiplicative reduction, this gives the right sign for semistable 
curves. The potentially multiplicative case has Wm{p) = wi(p)^™"'"^^/^. 

In the additive cases, the first author jj^j has used the work of Rohrlich [33| 
to compute the sign for p> 5. We get that"^" Wm{p) — u'i(p)'"^*p^/^, and wi{p) 
is listed in |2S1- For p = 2, 3 the value oiwi{p) is given^^ by Halberstadt |'23J, and 
our experiments for higher (odd) powers indicate that Wm{2) = 772^1(2)'^'" *^*J'-'/^ 
where 772 = —1 if V2{N) is odd and m = 3 (mod 8) and else ?y2 = +1, while the 
expected values of Wm{S) are given in Table 



Table 5. Experimental values for Wm{3) (periodic mod 12 in m.) 


$3 


1 3 5 7 9 11 


#3 


1 3 5 7 9 11 


$3 


1 3 5 7 9 11 


Cs,C4 

C2 


-F + -F -F + -H 
- + - + - + 




+ + + 

- + - -1- - -1- 


Ci XI Ci 
Qi XI C4 


-1- -1- - + 4- + 
----- -1- 



4.4 Computations 

From |5S], [51 Appendix B], or 16 , the assumption of the functional equation 
Am{E, s) = Wm,Am{E, m + 1 — s) allows us to compute (to a given precision) any 
value/derivative Am (E, s) in time proportional to Cm ~ V^m- Additionally, 
numerical tests on the functional equation arise naturally from the method. 

We follow § 7, p. 119ff|. Suppose we have Am{s) = WmA{m + 1 — s), and 
the dth derivative is the first one that is nonzero at s = k. Our main interest is 
in «; = u for m ~ 2u — 1 and k — v + 1 for m — 2v, and we note that d = for 
even m. Via Cauchy's residue theorem, for every real A > we have 



A 



id), 



d\ 



1 

2Td 



(•5) 



(-5) 



Am{z + k) dz 



yd+1 



The full conjecture^foUows from Langlands functoriality |28|. In the CM case, the 
functional equation follows from the factorisation @ and the work of Hecke \2,^ ■ 
Added in proof: A recent preprint 42 on Taylor's webpage shows the meromorphic 
continuation and functional equation for all symmetric powers for curves with j ^ Z. 
The work of Frohlich and Queyrut |22) and Deligne |13| might give a direct argument. 
Since we are assuming that m is odd, the exponent is just (m-|- 1)/2 unless <?p = C3. 
Note the third case in Table 1 of |23| needs a Condition speciale of C4 = 3 (mod 4). 



where S is small and positive and J, > is the integral along ^z = a. In the second 
integral we change variables z — > — 2; and apply the functional equation. Then we 
write K + X = m + 1, move both contours sufficiently far to right (say SRz — 2m) 
and expand Am in terms of the L-function to get 



A^i^l^ f + /,. , ,^^ V bmjn) 1 dz 

dl - 4„) ^" ^"^" + "^ ^^ n^^^ z<^^^ 2m A^ 

/ ,Nrf f ^r^\ / ^^•^bmin) 1 A^ dz 

J(2m) „_i 



n—1 



Thus we get that 
where 



^™ E ^^™(-; 3g;r) + (-1) V„c;^ E T^™(A; 



n— 1 n—1 



F-^r^.a-^- / 7m(^ + A^) dz 



(2m) z'^+^a;"^ 27ri' 

The FmifJ-', a;)-functions are "rapidly decreasing" inverse Mellin transforms. Note 
that we have lL (k) = /iL ('«)/7m(«;)C'm, and so can recover the L-value as 
desired. The parameter A allows us to test the functional equation; if we compute 
Am (n) to a given precision for A = 1 and A — 9/8, we expect disparate answers 
if we have the wrong Euler factors or sign Wm ■ 

We compute F^{fi;x) as a sum of residues at poles in the left half-plane, 
the first pole being at z = 0, following ^J. We need to calculate Laurent series 
expansions of the /^-factors about the poles. ^^ We let C(l) denote Euler's con- 
stant 7 « 0.577, and define Hi{n) — 1 for all n, and Hk{l) — J2i=i 1/* ^'^^ ^^^ ^i 
and recursively define Hk{n) = Hk-i{n) + Hk{n — l)/fc for n,k > 2. At a, pole 
z = — fc for k a nonnegative integer, we have the Laurent expansion 

r(z) . ±^^ (1 + ± H,{n){z + fc)") exp(f: ^^i^(. + kf 

and for k a negative integer (these only occur for a few cases) we can use the rela- 
tion zr{z) = r(z + 1) to shift. To expand r{z/2) around an odd integer z = — fc, 
we use the duphcation formula r{z) = r{z/2)r[^^)^^ to replace r{z/2) by 
a quotient of .T-factors that can each be expanded as above. The trick works in 
reverse to expand r(^^) about an even integer z = —k. We also have the series 
expansions for 2^ and 1/z about z — —k given by 



00 /i r\\T] -1 -, OC 



n— n—1 



^^ When Am,(Ei s) has a pole the factor 7m (s) has two additional linear factors (which 
are easily handled). But in this case it is better to use the factorisation J^. 



Since these F^{ii] x) functions are (except for CM) independent of the curve, 
we pre-computed a large mesh of values and derivatives of these functions, and 
then in our programme we compute via local power series. Thus, unlike the set- 
ting of Dokchitser |16| , we are not worried too much about the cost of computing 
F^(/x; x) for large x via a massively-cancelling series expansion, since we only do 
this in our pre-computations. For each implemented function we have its value 
and first 35 derivatives for all a; = i2'^/32 for 32 < i < 63 for k in some range, 
such as —3 < A: < 19. For sufficiently small x we just use the log-power-series ex- 
pansion. The choice of 35 derivatives combined with the maximal radius of a;/64 
for expansions about x implies that our maximal precision is around 35 x 6 = 210 
bits. When working to a lower precision, we need not sum so many terms in the 
local power series. Note that F^{ii;x) dies off roughly like exp(— a;-^/("'+^)), and 
thus it is difhcult to do high precision calculations for m large. 

To compute the meshes of inverse Mellin transforms described above, we used 
PARI/GP |2J|i which can compute to arbitrary precision. However, PARI/GP 
was too slow to use when actually computing the L-values; instead we used a 
C-based adaption of Bailey's quad-double package |2S], which provides up to 
212 bits of precision while remaining fairly fast.^"^ 

5 Results 

We tested the functional equation (via the above method of comparing the com- 
puted values for A = 1 and A = 9/8) for odd symmetric powers m = 2ii — 1 at 
the central point k = u, and for even symmetric powers m = 2v at the edge of 
the critical strip k = v + 1. We did this for all non-CM isogeny classes in Cre- 
mona's database [HI with conductor less than 130000; this took about 3 months 
on a cluster of 48 computers (each running at about 1 Ghz). 

We computed as many as 10* terms of the various L-series for each curve, 
which was always sufficient to check the functional equation of the third sym- 
metric power to about six decimal digits. ^^ In all cases, we found the expected 
functional equation to hold to the precision of our calculation. The results for the 
order of vanishing (at the central point) for odd powers appear in the left half 
of Table The right half lists how many tests^^ we did for other^^ symmetric 
powers (again to six digits of precision) . There are less data for higher symmetric 
powers due to our imposed limit of 10® terms in the L-series computations, but 
since the symmetric power conductors for curves with exotic inertia groups often 
do not grow so rapidly, we can still test quite high powers in some cases. 

^^ The SYMPOW package can be obtained from www.maths.bris.ac.uk/~inamjw 

^^ In about 0.3% of the cases, the computations for both the zeroth and first derivatives 
showed no discrepancy with A = 1 and A — 9/8; this coincidence is to be expected 
on probabilistic grounds, and for these cases we computed to higher precision to get 
an experimental confirmation of the sign of the functional equation. 

^^ We need not compute even powers when there is a lack of quadratic-twist-minimality. 

^^ We did not test the fourth symmetric power, as the work of Kim |2ti| proves the 
validity of the functional equation in this case. Since there is no critical value, a 
calculation would do little more than verify that our claimed Euler factors are correct. 



Table 6. Test-counts (right) and data for order of vanishing (non-CM isogeny classes) 



m 


Tested 


Order 


Order 1 


2 


3 


4 




m 


# tests 




m 


# tests 


1 


567735 


216912 


288128 


61787 


908 





6 


4953 


14 


26 


3 


567735 


262751 


287281 


16782 


905 


16 




8 


1259 




15 


1 even 


5 


46105 


22448 


23076 


569 


12 







10 


190 




15 


16 odd 


7 


3573 


1931 


1616 


25 


1 







12 


142 




16 


8 


9 


947 


542 


400 


5 










13 


5 even 




17 


3 odd 


11 


134 


51 


82 


1 










13 


30 odd 




18 


2 



Buhler, Schoen, and Top 6 already listed 2379bl and 31605bal as 2 examples 
of (suspected) 4th order zeros for the symmetric cube. We found 14 more, but 
no examples of 5th order zeros. For higher powers, we found examples of 3rd 
order zeros for the 5th and 7th powers, and 2nd order zeros for the 9th, 11th, 
and 13th powers, though as noted above, we cannot obtain as much data for 
higher powers. "'^^ We list the Cremona labels for the isogeny classes in Table [7| 

Table 7. Experimentally observed high order vanishings (non-CM isogeny classes) 



ord 



4th 
3rd 
2nd 



format is power:label(s) 



3:2379b 5423a 10336d 29862s 31605ba 37352d 46035a 48807b 55053a 

3:59885g 64728a 82215d 91827a 97448a 104160bm 115830a 

5:816b 23401 2432d 3776h 5248a 64801 7950w 8640bl 16698s 16848r 

5:18816n 57024du 7:176a 

7:128b 160a 192a 198b 200e 320b 360b 425a 576b 726g 756b 1440a 

7:15681 1600b 2304g 3267f 3600h 3600J 3600n 3888e 4225m 6272d 

7:11552r 158761 21168g 9:40a 96a 162b 324d 338b 11:162b 13:324c 



We also looked at extra vanishings of the 3rd symmetric power in a quadratic 
twist family. We took E as lla3:[0, — 1, 1,0, 0] and computed the twisted cen- 
tral value L^{Ed,'2.) or central derivative L'^{Ed,2) for fundamental discrimi- 
nants \d\ < 5000. We found 58 double zeros (to 9 digits) and one triple zero 
{d = 3720). A larger experiment (for \d\ < 10^) for 10 different CM curves found 
(proportionately) fewer double zeros and no triple zeros. 

Finally, we used higher-precision calculations to obtain the Bloch-Kato num- 
bers of equation Q for various symmetric powers of some non-CM curves of 
small conductor (see Table |Hll. More on the arithmetic significance of these quo- 
tients will appear elsewhere. In some cases, we were able to lessen the precision 
because it was known that a large power of a small prime divided the numerator. 

5.1 Other directions 

In this work, we looked at symmetric powers for weight 2 modular forms. De- 
launay has done some computations |llj for modular forms of higher weight; 
in that case, the work of Deligne again tells us where to expect critical values, 
and the experiments confirm that we do indeed get small-denominator rationals 
after proper normalisation. We looked at critical values at the edge and center of 

^'^ Given that we only computed the L- value of the 13th symmetric power for five curves 
of even sign, to find one that has a double-order zero is rather surprising. Higher- 
order zeros were checked to 12 digits; the smallest "nonzero" value was ~ 2.9 • 10~*. 



Table 8. Selected Bloch-Kato numbers for 





5th powers 


20a2 


2« 


37al 


2^ 


43al 


2^5 


44al 


2^'' 




Bth powers 


llaS 


2*5 


14a4 


2^3 


15a8 


210 


17a4 


212 


19a3 


2*3^5^ 


20a2 


2^73 


24a4 


2"/3 


26a3 


2^3 • 5 • 23 


26bl 


2^3 • 7=^ • 23 


30al 


f)15 oBrT 


33a2 


2^''3 ■5-7 


34al 


2"3359 


35a3 


2*3 • 7^31 


37al 


2^ 3'' 7 


37b3 2^3''467 | 


38a3 


2'^3*5 ■ 11 ■ 137 


38b 1 


2^3^5^13-31 


39al 


220 32 7 


40a3 


2^07/3 


42al 


2"3^7- 19 


43al 


2^3^1697 


44al 


22^5-31/3 


46al 


2^*5 ■ 23 ■ 30661 


50al 


2^5^73 


51al 


2^3^4517 



7th pow^ers 



24a4 
37al 
43al 



2^'' 7/3 
2133-5 
2^"^?, ■ 5 



lla3 
14a4 
15a8 
17a4 
19a3 
21a4 
24a4 
26a3 
26bl 
30al 
33a2 
34al 
35a3 
37b3 
38a3 
38b 1 
39al 
40a3 
42al 
44al 
45al 
46al 



9th powers 

212 

2"3*5 
2I6 

21^3^5 
2^325 
22«5 - 1 
2^^/ 9 



2I63555372 

2^*5- 107^167^ 
223355 . 72532 

2^53*5 



2" 5* 109^ 

2*»32 57* 



2^^^ - 223^241^ 



216319527.^32 
2"3i°53 14071^ 

l43 



50al 
54a3 
54b 1 



2«5 
2^3 ■ 5^ 



arious 


powers and curves 




10th powers 


lla3 


2"5 - 22453/3 


14a4 


2i'^335-6691 


15a8 


2^65 - 541 


17a4 


223327. ;^;l.227 


19a3 


2^3^47- 179-5023 


20a2 


24*53/3 


21a4 


22«3^5"29 


24a4 


2**^13/9 


26a3 


219357.47. 1787 


26bl 


219335273^^27-2102331 


40a3 


2^*5 - 683 


44al 


25^5 - 11 - 215447/3 


50al 


2"52«7/3 


52a2 


2**3^5 - 7 - 19 - 279751 


54b 1 


2143357 


56al 


2''«3^5- 11-71 


75c 1 


214528 . 31 . 41 . 61/9 


96bl 


2**197/3 


99al 


2IS331537. ;^3g7 




11th powers 


lla3 


226573 


14a4 


223365272-^^2 


15a8 


22911^2373 


17a4 


2263I1 


21a4 


23611^21173 


24a4 


2^^13745 


48a4 


2™ll73 


54b 1 


220341 


56al 


2^*1175 


72al 


268328592/5 



the critical strip, whereas we expect L-functions evaluated at other integers to 
take special values related to ii'-theory; see |HI3()I17I^ for examples. The pro- 
grammes written for this paper are readily modifiable to compute other special 
values. The main advantage that our methods have over those of Dokchitser 1^ 
is that we fixed the P-factors and the L-values of interest, which then allowed 
a large pre-computation for the inverse Mellin transforms; if we wanted (say) to 
compute zeros of L- functions (as with [SI), our method would not be as useful. 
Finally, the thesis of Booker takes another approach to some of the ques- 
tions we considered. The scope is much more broad, as it considers not only 
numerical tests of modularity, but also tests of GRH (§3.4), recovery of unknown 
Euler factors possibly using twists (§5.1), and also high symmetric powers (§7.2). 
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